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1 Introduction. 

Physics has often progressed very rapidly as the precision of measurements 
has increased. For instance, it was the precise measurements of Tycho Brahe 
which were instrumental in Kepler's deduction of the elliptic orbit, a result 
which later formed a cornerstone of Newton's universal law of gravitation. 
There is a minimum level of precision in the measurement of the planetary 
orbits, below which it becomes impossible to distinguish between elliptical or- 
bits with minuscule eccentricity, and circular orbits. Once data of the requisite 
precision to notice a difference was available, all that was left for someone to 
do is to put the pieces together, and voila you get Newton! 

Quantum Optics, plays an essential role in quantum metrology, a field in 
which the level of precision has increased exponentially over the past two 
decades. New techniques of increasing precision in quantum optics have in- 
creased the significant digits of some of the experimentally measured funda- 
mental constants by orders of magnitude. Quantum optics is so robust that it 



also of immediate use in testing theories of gravity and quantum field theory, 
for example, the L.I.G.O. collaboration are using a 2 km baseline Michelson 
Interferometer to search for gravitational radiation. 

These unheralded successes are currently pushing into new domains of ex- 
perimental precision, and we now have more direct access to the deeper layers 
of nature. Every week new quantum computing components are brought into 
being, by the shear effort of those working in the field. With each new switch, 
isolation mechanism, algorithm, etc., the goal of scalable robust quantum com- 
puting becomes more eminent. If we are successful in constructing quantum 
computers, the effect will be more revolutionary than anything before, includ- 
ing the classical computer and the internet. The vast expanse of Hilbert Space 
will then be in the throes of man. 

The fields of quantum optics and quantum computing are closely related 
to one another. Very often breakthroughs in quantum optics are implemented 
in quantum information processing, storage and quantum communication de- 
vices. For example, two ways in which cavity QED techniques may be used to 
to perform quantum computations are (from [8]) 

1. Quantum information can be represented by photon states, with atoms 
trapped in cavities providing the non-linear interactions between pho- 
tons, necessary for entanglement. 

2. Quantum information can be represented by atoms in different states, 
where photons are used to communicate between the different atoms/states. 



Any realization of these schemes would at some point have to address the prob- 
lem of precision control of population transfer, as a means to generate single 
photons. Such precision is a per-requisite for realizing any completely quan- 
tum technology, that is, any technology based on computational components 
whose functionality depends a'priori on quantum non-linearities, an example 
of which is entanglement. 



2 Field Quantization. 

This treatment of the field quantization will closely (but not exactly) follow 
chapter 2 of Gerry et.al., given in [1]. In order that we understand the inter- 
action of quantized modes of the electromagnetic field with "atoms", (whose 
definition will, for the moment remain general; we will define an atom to be 
any bound sate of electrons in a potential V (r).) we must first understand the 
properties of the quantized fields themselves. In the following we begin with 
the simple case of a single mode field confined to a 1-d cavity. This clearly 
represents an idealized situation, but we will later generalize to the case of a 
multimode field in some three dimensional cavity. 

2.1 Single mode field. 

We begin as always, with the one-dimensional square well, but in the con- 
text of quantized modes of the electromagnetic field, which will be relevant for 



our later analysis of quantized modes of optical cavities , etc. One fruitful and 
interesting scenario to investigate for our purposes, is the case of a radiation 
field confined to a one dimensional cavity free of sources(i.e. there are no cur- 
rents, charges, or any dielectric media in the cavity), oriented along what we 
choose to be the z axis, with perfectly conducting walls at z = and z — I, 
therefore the transverse electric field must vanish at the boundary. 

Recall that in SI units, the source-free Maxwell equations, which our single 
mode field must satisfy, are 

VxE = ^- (2.1) 

at 

<9B 

V x B = fi e ^- (2.2) 
at 

VE = (2.3) 

V • B = (2.4) 

We will assume that the field is polarized in the x-direction i.e. E (r, t) = 
e x E x (z, t) and hence B (r, t) = e y B y (z, t). If we identify q (t) as the canonical 
position as denned in the Hamiltonian formalism and similarly identify q (t) 
as the canonical momentum, then the solution for the components is 



E x (z,t) = [ — ) q(t) sin (kz) (2.5) 



and 

,2 \ 1/2 



= (^r) (^r) 9 (*)«»(**) (2.6) 



where the wave-number k is related to the frequency ou by k — ui/c. Moreover, 
the boundary conditions on the electric field at the interface of the perfect 
conductor at z = and z = L, constrain the values of k to be 



/rar. 

/■• (— ) ■ " = 1-2. 



(2.7) 



and therefore the allowed frequencies are 



/mr\ 

w = c ^— J , n = 1,2, 



(2.8) 



We can invert the 2 equations giving E x and B y in terms of the canonical 
position and momentum q (t) and q (t) , and obtain the expressions for q (t) 
and q (t) in terms of E x and B y , namely 



q (t) = E x (z, t) [ J esc (kz) 



(2.9) 



q{t) = B y (z,t) 



sec (fc^) 



(2.10) 



From these expressions, it is apparent that the Hamiltonian for the field is 



H 



dV 



e E 2 (r,t) + — B 2 (r,t) 



(2.11) 



Now, 



e E 2 (r, t) = e E (r, t) ■ E (r, t) = e (e x ■ e x ) E 2 X (z, t) 



= e E 2 x (z, t) 



and similarly, 



Therefore 



From (1) we have 



and 



Therefore (6) becomes 



1b 2 M) = -V 2 M) 

Ho Ho 



H=\fdV 



e E 2 x (z,t) + -B 2 y (z,t) 



2co 2 

e E 2 x (z,t) = —q 2 (t) sin 2 (kz) 



—B 2 (z,t) = ^p 2 (t) cos 2 (kz) 
Ho v 



H= i rdv ^ 2 ^ sin2 ^ + p2 ^ cog2 ^ 

2 J V 



Since, 



and, 



1 + cos 2x 



cos x = 



1 — cos 2x 



sin x = 



we may write the Hamiltonian as 



H= - I ^- [cu 2 q 2 (t) (1 + cos2fcz) + p 2 (t) (1 - cos 2kz)] 
2 J V 



Now the cosine terms drop out of because of the periodic boundary conditions 
and therefore, 

= I(p 2 +a;V) (2-21) 

and so the system is equivalent to harmonic oscillator with unit mass, (q (t) = 

P(t)). 

Now that we have the canonical momentum and canonical position, is is 
relatively easy to quantize the field by replacing the variables H with H and 
q (t) and q (t) with the hermitean (observable) operators q and p, respectively. 
Moreover, we must require that the observables obey the canonical commuta- 
tion relation 

[q,p]=ihl nxn (2.22) 
which will write from here on out simply as44 



[q,p]=ih (2.23) 

with the n x n matrix identity operator I nxn implied. Having promoted q and p 
to operators, we are thereby led to the operators for the electric and magnetic 
fields 

2c^ 1/2 



E x =[ — ) q sin (kz) (2.24) 



and naturally, the Hamiltonian operator becomes 



H = \ {p 2 + co 2 ?) (2.26) 



Now we define the non-hermitean creation, a), and annihilation, a, opera- 
tors as follows[l]: 

Vm^a^ = (coq - ip) (2.27) 



Denning 



and 



V2frwa = (uq + ip) (2.28) 



S = (hu/Ve ) 1/2 (2.29) 



B = (fi /k)(e hLo 3 /V) 1/2 (2.30) 



it follows that we can write the operators for the electric and magnetic fields 
as [1]: 

E x (z, t) = £ (a f (t) + a (t)) sin (kz) (2.31) 

B y {z, t) = iBo (a f (t) - a (t)) cos (kz) (2.32) 

From now on we will suppress hats, ' , on operators and just write a) = a\ 
a = a, E x = E x , etc. The benefit of working with creation and annihilation 



operators is that we are allowed to utilize the simplicity of their algebra. 



[ a , a) ] = aa) — a) a = 1 



(2.33) 



These commutation relations allow us to write the Hamiltonian operator as 
[2] 

H = fku (at a + (2.34) 

In the Heisenberg representation, a general operator O will obey Heisenberg's 
equation of motion [2] 



dO _dO i 
dt dt h 



H,0 



(2.35) 



Which in the case that O does not depend explicitly on the time coordinate, 
becomes 



dO__i_ 
dt h 



H , O 



(2.36) 



Therefore for the creation and annihilation operators we have the following 
time evolution equations[l] 

An} 

(2.37) 



dt 



= iooa 



da 
~dt 



= —iLoa 



(2.38) 



which implies that 



a t ( t ) = a t (o) e ^ 



(2.39) 



and 

a(t) = a (0) e~ iu}t (2.40) 

We may expand e~ luJt as 

e- wt = l-iut--^- + i-y- + ... (2.41) 
which allows us to write a (t) as 

(22 3 3 \ 

A useful combination of operators will be a+a = n , a combination known as the 
number operator. If applied to the n th eigenstate of the Hamiltonian \n) (we 
will later come to identify \n) as the n photon state), 

(J a \n) = n \n) (2.43) 

this operator gives the value n of the eigenstate occupied. The energy eigen- 
value problem can then be written as 



H ;/) = h,s -in ■ ^ ) \n) (2.44) 



[n+-)\n) = E u //) (2.45) 



E n = Hco I n + \ ) . (2.46) 



Therefore 

E„ = fiuj I n 4- 

2 

Where E'o is the ground state energy, since is the lowest value which may be 
taken by n as can be seen from acting on the state |0) with the annihilation 
operator a. 

a|0) = (2.47) 

Since 

a\n) = y/n\n - \) . (2.48) 

and 

a) \n) = Vn + 1 |n + 1) . (2.49) 

It follows that any arbitrary eigenstate \n) can be written in terms of the vac- 
uum state as (e.g. [1]) 

\n) = ±jL- 10) . (2.50) 
vn! 



The states \n) form a complete basis for the Hamiltonian H, and are orthonor- 
mal 

(m\ \n) = 5 mn (2.51) 

The non vanishing matrix elements of the creation and annihilation operators 
are 

{n-l\a\n) = y/n (2.52) 



(n + 1| a) \n) = y/n + 1 



(2.53) 



2.2 Quantum fluctuations of the single mode field. 

Recall the operator for the electric field given in (8), 

E x (z, t) = S (a 1 " (t) + a (t)) sin (kz) . (2.54) 

The eigenstates of the Hamiltonian \n), do not form a basis for the operator 
E x (z,t). This is implied by the fact that the number operator n = a) a, which 
does commute with the Hamiltonian, does not commute with the electric field 
operator E x . 

To see this, let us first calculate the average field (E x ) 

(E x ) = (n\ E x (z, t) \n) = S [(n\ a \n) + (n\ a t \n)] sin (kz) (2.55) 
which, by (10) and (11), become 

= €o[0 + 0]sm(kz) = (E x } = (2.56) 

that is, the average field is zero. 

The energy density of the field E = e x E x is proportional to the mean square 
of£ x ,[l] 

(E 2 X ) = 2£l sin 2 (kz) (n + ^ (2.57) 



The variance is defined as [1] 



<(A^) 2 ) = {El) - {E x ) (2.58) 

i.e, its the mean square of the standard deviation, which for the eigenstate 
\tp n ) becomes 



AA",. - \/ (/:;} ( E x ) = J25 2 sin 2 (A;*) U + ^ ) (2.59) 



that is, 

/ 1\V2 

AE X = V2£ Q sin (kz) ( n + - J (2.60) 
It is interesting to note that even for n = we have 



A.E,. = v^o sin (kz) = £ sin (Jfez) (2.61) 

V2 

these are called the vacuum fluctuations of the field, since they correspond to 
the eigenstate of the vacuum |0), the state with zero photons. [1] 

In the case of the electromagnetic field confined to a ID cavity, the eigen- 
states of the Hamiltonian, namely|n), correspond to states of photon number n 
. One important fact is that the number operator n = a^a and the electric field 
operator E x do not commute, 

[n, E x ] =S sin (kz) (a f -a) (2.62) 



The generalized uncertainty relations state that for any two operators A and 
B satisfying [ A , B ] = C , it follows that the product of the uncertainties of A 
with that of B obey the inequality 




(2.63) 



It follows therefore, that the number operator and the electric field obey the 
following uncertainty relations 



where < < 2n is the phase angle associated with the creation and anni- 
hilation operators. In fact, it will turn out, that the situation is not actually 
quite that simple. It turns out to be a very slippery task to define a unique 
phase operator, and in fact is not possible in general [5] ,[6], [7]. It can be 
shown, however, that for proper definitions of the phase,namely those given 
in [4], that the photon number states \(p n ) have a uniform phase distribution 
(A0/An) « constant for < < 2n. For more on the number phase uncertainty 
relations see [3], [4] and the references given in [1] . 




(2.64) 



This implies a number-phase uncertainty relation [1] 



An A0 > 1 



(2.65) 



2.3 Multimode fields. 



In free space in the absence of any sources, the source free Maxwell equa- 
tions are still valid, (joking but obviously true). We write the electric and mag- 
netic fields in terms of the vector potential A (r, t) which satisfies the wave 
equation [1] 

1 <9 2 A 

v2A -4^r = ° < 2 - 66 > 

& at 1 

and we choose the Coulomb gauge condition( which will become useful later on) 

V • A = (2.67) 
The electric field is then given by[l] 

E = ~ (2.68) 

at 



and the magnetic field is 



B = V x A (2.69) 



As long as I > | , we can model free space as cubic cavity, with sides of 
length L, therefore we may impose periodic boundary conditions on the faces 
of the cube[l]. This allows us to deal with the mathematically simpler case 
of having a denumerably infinite set of normal modes, rather than a non- 
denumerably infinite set of modes [1]. We require plane waves in the Xi di- 



rection, where (i = 1, 2, 3) and (xi = x, x 2 = y, x 3 — z), to satisfy the condition 



e ik Xi Xi _ e ik Xi (xi+L) ^ 70) 

which leads the following conditions for the direction numbers k Xi 

(2.71) 



2vr^ 




~Lj 


) m x 


2tt^ 






) m v 


27^^ 




T > 


) m z 



k y = m y (2.72) 

/9tt\ 

(2.73) 



where 



Now, the wave vector 



m x = m y = m z = 0,±l,±2,... (2.74) 



2tt 

k = (k x , k y , k z ) = — (m x , m y , m z ) (2.75) 



Moreover, k = ||k|| = Vk • k = uok/c . Distinct normal modes of the fields are 
specified by distinct sets of integers (m x , m y , m z ). Therefore, the total numbers 
of modes in the interval (Am x , Am y , Am z ) is [1] 



(-' 



Am = Am x Am y Am z = 2 — Ak x Ak y Ak z (2.76) 



taking into account a factor of 2 for the two independent polarizations. In the 



limit that L — > oo, Am — > dm and we have (V = L 3 ) 



dm = dm x dm y dm z = ( ^— ^ ) dk x dk y dk z (2.77) 



going to spherical coordinates this is 

k = (k x , k y , k z ) = k (sin 9 cos <j), sin 9 sin 0, cos 9) (2.78) 



and therefore 

dm = I 

or using k = Uk/c we can write this as 



dm= ( ^-)k 2 dkdn (2.79) 



dm = ( ^du k dn (2.80) 



Integrating over the solid angle Q, we obtain 

V 



dm = —k 2 dk = Vp k dk (2.81) 
ir 2 



Where pk = k 2 /n 2 . We may also write this for duk as 



<jj 2 

dm = V^rduk = Vp M duj k (2.82) 



for which 

p (cu k ) = ul/ (ttV) (2.83) 



Having "pasted" the cubic grid on our space, we may proceed to expand the 
vector potential as 

A (r, t) = J2 e ks [A ks (t) e* kr + A* ks (t) e"^] (2.84) 

k,s 

where A ks e C is the amplitude of the field and e ks elisa polarization vector 
[1]. Moreover the sum over k is the sum over the distinct sets of integers 
(m x , m y , m z ), and the sum over s is the sum over the two polarization directions 
[1], which must obey the orthonormality relations 

e ks ■ e ks / = Skk'Sss' (2.85) 

The Coulomb gauge condition requires that k ■ e ks = 0, which is known as the 
transversality condition [1]. 

The wave equation and the Coulomb gauge lead to the following relations 
for the amplitudes A ks : 

^ + oolA^ = (2.86) 
at z 

The solution to this differential equation is 



A ks (t) = A ks e 



(2.87) 



(A^s (0) = A ks ). Thus, the electric and magnetic fields become [1] 



E(r,t) = i^ cu k e ks [A ks e^ r -^ + A* ks (t) e -*0"-«*')] 

k,s 



B (r,t) = J J> fe x e ks ) [A ks e^ r -^ + A* ks (t) e -*( k -«**)] 



k,s 



The energy of the field is 



H 



dV 



e E ■ E + — B • B 



Now [1], 



x e k , 



k 

Iki 



x e ks / = 5 SS > 



and 



k 

in 



k 

IkT 



Taking our periodic boundary conditions into account, we have 



/ 



e ±i(k-k')- w = 5kk/V 



Therefore, the contribution to H from the electric field is 



\ J dV [s E ■ E] = e V "l^sAl s - R 

k,s 



The contribution from the magnetic field is 

I f dV ± B ■B = e Vy2 ^ s A* ks + R (2.95) 
1 J , 



where, 



R = \e a V "fa* ■ e^ s , [A ks (t) A^ s , (t) + A* s (t) A\ s , (t)} (2.96) 



2 



Therefore, the total energy in the field is 

H = 2e V (t) A* ks (t) (2.97) 

k,s 

but since A ks (t) = A ks e~ luJkt , this may be written as 

H = 2eoVj2^tA^A* ks (2.98) 



k,s 



2.4 Quantization of the multimode field. 

We may quantize the field by introducing the canonical position and mo- 
mentum operators g ks and p ks , respectively through the definitions 



Iks = - — K= [uJkqks + ipks] (2.99) 
2u k ^e V 



A L = -z — T=f7 \- Wk< i^ ~ (2.100) 



in which case the Hamiltonian becomes 



(2.101) 



as it should. 

The canonical variables obey the canonical commutation relations 



[<?k s , <7kv] = = [pks,Pkv] 



(2.102) 



(2.103) 



Just as we did for the single mode field, we may define the creation and anni- 
hilation operators for the multimode fields 



(2.104) 



y/2ftw k al s = uj k qks - ipks 



(2.105) 



The creation and annihilation operators obey the following commutation 
relations 



[Oks , Ok'. 



fl L ) a t's' 







(2.106) 



Oks j Ok' 



8kk>8 ss >5 ( k - k 



(2.107) 



Just as was the case for the single mode field, the number operator for the 



mode ks is n ks = a ks a ks , and the Hamiltonian is 

# = («L«k, + ^) (2.108) 

= ^^ fc (n ks + ij (2.109) 

Each mode is independent of all the rest and has the eigenstates |n ks ). If 
we let j denote the j th mode kjSj, then we may write the n th photon number 
state of the j th mode as 



r'" 



■j 



\M) = U^rr rf'O (2.H0) 

3 v J 

The energy eigenvalue equation is then 

H\{n J }) = E\{n J }) (2.111) 

where [1] 

E = J2^U + l) (2.H2) 

j 

A multimode photon state is the tensor product of all of the individual mode 
number states, that is 

|ni,n 2 ,n 3 ...) = \ni) <g> \n 2 ) <8> |n 3 ) <g) ... (2.113) 



(2.114) 



The number states are orthogonal, that is 



3,3 



The multimode vacuum state is 

|{0j}) = |0) <g> |0) <g> |0) <g> ... (2.116) 

The action of the creation and annihilation operators on j th mode of the multi- 
mode photon number state are given by 

a S \{n s }) = y /fq\{(n-l) j }) (2.117) 



a] I {n 3 } ) = v^TTT | { {n + 1 ) . } ) (2.118) 
Quantization requires that the amplitudes A ks become the operators: 



1/2 



<*• (2 ' 119> 



Which therefore allows us to define a vector potential operator as well as elec- 
tric and magnetic field operators, which are, respectively, 



A(r,t) = (^Tv) ' Gks [ ak ^ i(k ' r " U ' fc * ) + aLe" <(k-r " Wfct) (2-120) 

k,s 



E(r,i) 



k,s 



2e n V . 



1/2 



a ks e 



i(kT-wfct) 



i(k-r-aj fe t) 



(2.121) 



B(r,t) = ±J> (^ X M (k ^ - aLe- i(k - r - Wfct) ] (2-122) 

k,s VII / V / 

where the operators ciks = Ok s (0) form a basis of the Heisenberg representation. 
The time dependent creation and annihilation operators are, respectively 

«LW = «L(0)e^ (2.123) 

a ks (*) = a ks (0)e-^ (2.124) 

It can be seen that the magnetic field is weaker than the electric field by a fac- 
tor of 1/c, which is what we should expect and is a reassuring sign we haven't 
gone off track. The magnetic field couples to the spin magnetic moment of the 
electrons which is negligible for the aspects of quantum optics which we will 
investigate. 

Some interesting things to note are: 

• A single mode plane wave has electric field components given by 

/ nil \ 1/2 

E (r, t) = if e x lae^'-^ - a) e'^^} (2.125) 

\2e VJ 1 J 



Quantum optics often works in the domain of optical radiation, whose 



wavelength A is on the order KrA.Therefore in these situations we may 
approximate 

e ±ik - r « 1 ± ik ■ r (2.126) 
since, in these situations it is true that 

A = _L > |r atom | (2.127) 
2n |k| 

and hence the electric field can be expanded as 

/ hi \ 1/2 

E (r, t) « E (t) = i [ae~ wt - ate™] (2.128) 
which is known as the dipole approximation . 



3 Thermal Modes. 

Consider a single mode field in thermodynamic equilibrium with the walls 
of a cavity of absolute temperature T. The density operator p for the system is 
[2] 

p = Z- l e- H ' kT (3.1) 

where H is the Hamiltonian, k is Boltzmann's constant and Z is called the 
partition function, and it is introduced as a normalization factor, in order that 



trace of p be one[2]. That is 



Z = Tr {e- H ' kT } . (3.2) 

The density matrix for the system is diagonal in the Hamiltonian eigenbasis 
\(p n ) since it is in thermodynamic equilibrium. The diagonal matrix compo- 
nents, gives the population of the stationary state | </?„), which in this case are 
all the same. They are: 

Pth = Pun = (n\ Z-\- u l kT \n) (3.3) 

= Z- 1 e- J5 »/ fcT . (3.4) 
Since E n = huo (n + |), the partition function becomes[l] 

Z = exp {-hu/2kT) ^ exp {-hun/2kT) (3.5) 

n 

Since , exp (—ftuj/kT) < 1 we may sum the series 

£ exp (-hu/k B T) = 1 _ exp( _ 1 _^ /A;T) (3-6) 

Therefore 



exp (—hu/kT) 



(3.7) 



1 — exp ( hjj/kT) 

The off-diagonal terms vanish, hence there are no coherences between sta 



tionary states, and 



Prim 



(n 



Z- l e- H l kT \m) 



(3.8) 



z -l e -E m /kT ^ = q 



(3.9) 



We observe that in thermodynamic equilibrium, the populations of the sta- 
tionary states decrease exponentially with the energy. Since there are no co- 
herences, the system in this case may be considered to be a statistical mixture 
of the states \n). 

The probability that the thermal mode is in the n th thermally excited state 

is 



P n = (n\pth \n) 



(3.10) 



exp(-E n /kT) 



(3.11) 



£ n exp(-£ n //cT) 



The density operator may be written as[l] 



'OO 




(3.12) 



n'=0 n=0 




(3.13) 



n=0 



OO' 




(3.14) 



n=0 



The average photon number of the thermal field is [1] 

n = Tr(n ^ ) = exp(^/A;T)-l (3 ' 15) 
from which it follows that for kT > Hu — >■ 

n « fcT/Tju; (3.16) 

While for hu ^> kT ^ 

n « ftw/A;T (3.17) 

4 The interaction of atoms and electromagnetic 
waves 

The Hamiltonian for a system consisting an electron bound to an atom in 
the presence of external fields is 



H (r, t) = — [P + eA (r, i)] 2 - e$ (r, t) + V (r) (4.1) 
2m 

The gauge invariant electric and magnetic fields are given by 

E (r, t) = — V$ (r, t) — — (4.2) 

B(r,t) = V x A(r,t) (4.3) 



Gauge invariance means that these fields are invariant under the gauge trans- 
formations 

$'(r,t) = $(r,t)-^^ (4.4) 
A'(M)=A(M) + V X (M) (4.5) 

Therefore the equation governing the time evolution of the system in the 
Schrodinger representation, is the Schrodinger equation 

(r t) 

H(T,t)*(i,t)=ih ^ ' . (4.6) 

In quantum mechanics all operators are invariant under obey a global U (1) 
similarity transformation; for some operator A, A' = UAU\ where U is some 
unitary operator. Moreover, all state vectors are invariant under multiplica- 
tion by a common U. This essentially means that given some quantum me- 
chanical representation of system with all of its operators, and states, etc., 
we may obtain an equivalent description of that same system if we simulta- 
neously transform all of the states and operators of the theory in the manner 
prescribed above. The resulting transformed theory will lead to all of the same 
results as the original theory. The usefulness of this fact is that a particular 
operator may take on a more tractable form in the transformed theory. There- 
fore, we may exploit this fact to simplify the Hamiltonian. It will prove useful 
for us to define the unitary operator R which takes us to another representa- 



tion & (r, t) of the eigenstate * (r, £), by the action of R on ^ (r, t), namely 



m (r, i) = (r, t) . (4.7) 



The transformed Hamiltonian obeys its own Schrodinger equation 



H' (r,t)V (i,t)=ih 9<if £' t) (4.8) 



where [1] 



Choosing 



H' (r, i) = + ih^rrf (4.9) 



i? = exp (-iex (r, t) /ft) (4.10) 
which amounts to choosing the Coulomb gauge, we have 



2m 



2 



P + eA(r,t) - e$ (r, i) + V (r) (4.11) 



It is important to note that we will be working in the Coulomb gauge, which 
is not relativistically covariant, but for which $ (r, t) — and V • A = ( 
the transversality condition), therefore the radiation field is completely deter- 
mined by the vector potential . If there are no sources near the atom, then A 
satisfies the homogeneous wave equation 



(4.12) 



whose solution has the form 



A = Aoe i(k - r - Wfct) + Aje- i(k - r - Wfct) (4.13) 

In the Coulomb gauge the radiation field is completely determined by the vec- 
tor potential as can be seen from the Hamiltonian 

H (r, t) = ^ + — A • P + —A 2 + V (r) (4.14) 

The transformed Hamiltonian becomes[l] 



H' = ^-[P + e(A + V X )] 2 - + V (r) (4.15) 



Since |k| = 27r/A, it follows that for |r|~ 5a ( Bohr radius) and A ~ 500nm 
(optical radiation) , k • r l,Thus we may invoke the dipole approximation 
which gives the first order interactions and which also implies that on length 
scales ~ |r| atom , the vector potential is locally uniform, and we may make use 
of the fact that A (r, t) ~ A (t). If we choose our gauge function to be 

X (r,t) = -A-r (4.16) 



it follows that 



Vx (M) = 



-A(t) 



(4.17) 



which means that 



p2 

H'(r,t) = — + V(r) + er-~E(t) (4.19) 
2m 



we recognize the quantity d = — er as the electric dipole moment, and we may 
write 

H' (r,t) = H -d-E (t) (4.20) 

4.1 Interaction of an atom with a classical dipole field 

Let us begin with the case of a classical field of frequency u, given by [1] 

E (t) = E cos (ut) (t) (4.21) 

where 

e (0 = { (4.22) 




just means that the field is turned on at a time t = 0. We can study the in- 
teraction of an atom with this field by using perturbation theory. Expanding 
to first order just amounts to using the dipole approximation, k • r 1 , which 
we have seen previously is satisfied in the case of atoms interacting with a 
classical electromagnetic field. 

Given an atom, in some initial state we can expand of the atomic state 
of the atom for all t > 0,in a basis of uncoupled atomic states \k) , which span 



the space of H mt , 

\^(t)} = j2 c k( t ) e ~ tEkt/h \ k ) ( 4 - 23 > 

k 

where the amplitudes C k (t) are normalized such that 

£l<7*(*)| 2 = l- (4.24) 

k 

Now, working in the Schrodinger picture, the atomic state at time t , \ip (t)), 
must obey the time dependent Schrodinger equation, which is: 

ih^^ = (H + ir a )\il;(t)) (4.25) 

where, in the dipole approximation, as we know from the last section, H mt = 
— d-E (t). Substituting our expression for \ip (£)) into the Schrodinger equation, 
and then multiplying from the left by (/| e~ tEkt / R ,we ( denoting time derivatives 
with dots ' ) obtain a set of coupled first order differential equations for the 
amplitudes 

Ck(t) = - l -Y,C k (t) {l\H mt \k)e^ 1 
n k 

Where, 

u lk = (E l - E k ) /h, (4.26) 

are the transition frequencies between atomic states |/} and \k) . In order that 
we may solve these equations we must also subject them to the condition, that 
the initial atomic state is \i) , which implies that C k (0) = 1. As the state evolves 



in time, the initial state will transfer to other other |/) with a probability given 
by 

P^f (t) = \C,(t) | 2 . (4.27) 

To further simplify the task of solving this set of coupled equations analyti- 
cally we expand the amplitudes as a power series in some coupling parameter 
< A < l(which measures the strength of the interaction relative to scale at 
which our theory breaks down, A = 1). 

Q (t) = cf } (t) + ACf } (t) + A 2 cf } (t) + ... (4.28) 

Inserting the expression for C\ (t) into (IV. 1) we obtain a recursion for the 
n th amplitude 

c i n) (0 = -jX ^ {t) (/| H * lk) ^ (4 ' 29) 

k 

which leads to a coupled set of equations for all of the C" (t), which up to second 
order are given by 

cf (t) = (4.30) 

w = 4 £ ^ 0) w <*i ^ i*> e<wttt (4 - 3l) 

cf > w = 4 E c k ] (*) ci ^ e<w,fct (4 - 32) 

The only surviving terms in the sum are those for k — % . Therefore the first 



order amplitude becomes, upon integrating on time 



Cf (f) 



if'""' 



int e iu, fi t> C (0) ^ 



(4.33) 



using the recursion relation, inserting this value for (t) into the equation 
for (jj® (t) and integrating on time enables us to find (t), (see [1]) 



(4.34) 



[ dt' I dt"Hf (*') e^'Hlf (t") e^'cf (f) (4.35) 



o ./o 



The total transition probability as a function of time, for a transition from 
state \i) to a state |/) is: 



p^ f (t) = cf it) + c) L> it) + cf (t) + ... 



(i) 



y(2) 



(4.36) 



We have up to now, neglected taking account of the fact that the dipole mo- 
ment operator d only has non-vanishing matrix elements for states of opposite 
parity. Taking this into account we see that up to first order cf (f) = 1, so 
that 

• nt 

Cf\t) = ~\] dt'H'/fe^' 

' e K"+"f)t _ A ( e -i(u+u f )t _ { 



(4.37) 



1 

2h 



d " E ) fi 



(U + UJfi 



{OJ - UJfi) 



(4.38) 



When the radiation frequency u is near the atomic transition frequency w/j, we 



will have resonance, and in this case we may neglect the first "anti-resonant" 
term since the second term will clearly dominate. This is the so called rotating 
wave approximation. With this in mind, the first order transition probability 
becomes [1] 



where we have introduced the notation A = u — u/i, which is known as the 
detuning parameter. 

4.2 Interaction of an atom with a quantized dipole field 

Earlier, while working in Heisenberg representation, we found that the 
quanta of a single mode electric field, in the absence of sources of any kind, 
were given by 




(4.39) 




(4.40) 



switching to the Schrodinger representation, this becomes 




(4.41) 



The free Hamiltonian is 



H = H atom + H f 



(4.42) 



where 

H atom = p2/2m + V (r) (4.43) 

and 

H fteld = hwcta (4.44) 

are the source free Hamiltonians of the atomic system and the field, respec- 
tively. We have suppressed the vacuum energy in our expression for H f%eld 
because it does not contribute to the dynamics. The interaction Hamiltonian 
is 

H M — — d -E(t) — -i {^—) (d ■ e) (a - t) (4.45) 
= d • S (a f - a) (4.46) 

where, 

S = i{hu/2e V) 1/2 e (4.47) 

We have thus quantized the atomic system as well as the field system. 
If we wish to combine the distinct atom and field systems into one, atom- 
field system, we must remember that the state space of the atom-field system 
will in general be a linear superposition of the eigenstates of H atom and Hf ieid . 
Consider an atomic system in the initial state \a) . If we combine this atomic 
system with the field system which initially contains n photons, then we will 
have the atom-fields system which is initially in the state 

\i) = \a) \n) . (4.48) 



Since the interaction Hamiltonian H mt is proportional to (a> — a), it follows 
that for the n th eigenstate \n), the only non-vanishing matrix elements of H mt 
(in the atom-field eigenbasis) are the following 

(H mt ) = (fi\ Hint \i) = (d • £ ) ba (b , m\ (a* - a) \a , n) = (4.49) 

i=l,2 

= (d • £ ) ba (y/n + l<J„,.„+i - a/^W-i) • (4.50) 

where 

(d^o) 6a = (&|d|a)^ (4.51) 

The quantity (6| d \a) = d ha gives the transition dipole moments between the 
states 1 6) and | a) .Therefore the interaction Hamiltonian couples the n th state 
to either the n+ 1 or n — 1 state. In fact H mt induces a transition from the initial 
state of the atom-field system \i) to the state = \b) \n — 1) by absorption of 
a photon or to the state |/ 2 ) = \b) \n + 1), by the emission of a photon. The 
energies of these states are [1] 

\i) = | a) |n) -H- = E a + n/ja; (4.52) 

\h) = \b) \n - 1) % = E 6 + (n - 1) ^ (4.53) 

|/ 2 > = |6>|n + l) ^ % = E b + (n + 1) ^ (4.54) 

where E a and E b are the energy eigenvalues of the respective atomic states \a) 
and |6). 



Let us compare the results of the semi-classical versus the quantum treat- 
ment of this problem. In both cases, absorption is forbidden in any state for 
which n = (zero photons in the system), for obvious reasons. However, for the 
quantum case of emission, even if n = transitions may occur known as spon- 
taneous emission a phenomenon with no semi-classical analog. In cases where 
n > 0, we then speak of the stimulated emission of an additional photon. In 
the classical case, if you start with no field, i.e no photon, then you will never 
have a photon later, but a photon later is almost a certainty in the quantum 
case. The matrix elements of the interaction are in the case of adsorption [1] 

(/^ H int \i) = (b,n-l\ H int \a,n) (4.55) 

= -(d-£ ) ba ^ (4.56) 
and for the case of emission, are 

(f 2 \H int \i) = (b,n+l\H mt \a,n) (4.57) 

= (d-£ ) ba V^TT (4.58) 

where, just as before, 

(d-£ ) ba = (b\d\a)-E < 4 - 59 ) 

Fermi's golden rule tells us that the rates of emission and absorption are 
proportional to square modulus of the matrix element coupling initial \i) and 
final states , |/ 2 ) , which in the case of a single mode (monochromatic) field 



coupled to an atom, whose final state space is spanned by \a), \b), The transi- 
tion matrix elements are given by: (see [1]) 



7T 

[/] 



W^ m = | £ 1 ^^-S (« - . (4.60) 



Moreover, since, 

(d-£ ) ba = (b\d\a)-S (4.61) 



this becomes, 



W 71 \( b \ d \ a ) '^0| 2 X/ x ^ CO\ 

wi-i/i = ^ 5 ^ ~ w ^ ■ (4 - 62) 

[/] 



Therefore the ratio of these rates is given by 



\(h 


jjint 


*)| 2 


\(h 




°>l 2 



n 



(4.63) 



For the time being, let's just focus on two atomic states \a) and \b), and 
ignore the rest. This allows us to write the write the state vector as 

(t)) = Q (t) \a) \n) e -iE*t/h e -ir** + Cfi ^ ^ | n _ ^ e -iB 6 t/R e -i(n-l)«t 

+ C /2 (t) 1 6) |n + 1) e - lEbt/h e- l{n+1)u)t (4.64) 

since we already said that the initial state was \ip (t)) = \a) \n), so therefore we 
have: 

Q (0) = 1 (4.65) 



C/i (0) = C f2 (0) = (4.66) 

We can use perturbation theory just before to obtain the probability ampli- 
tudes for all times, t > . 

C}? (t) = ~ dt' (fi\H mt \i) e^n^'cf (t') (absorption) (4.67) 
" Jo 

C$ (*) = --/ dt' (f 2 \ H int \i) e^*'cf } (t') (emission) (4.68) 
' " Jo 

where Uf U = (E fl — Ei) /h and u f2i = {Ef 2 — E^) jh . Therefore the probability 
that the atom under goes a transition to the final state \b) is the sum cj 1 ^ = 



Cfi + , i.e. 



Cf ] == \ (d • E ) a6 



(4.69) 



where oob a = (E b — E a ) /h. If the initial state \a) happens to be the excited state 
then bj ha < 0. In this case, for radiation frequencies u ~ (— u ba ) , and we get 
spontaneous emission. 

4.3 The Jaynes-Cummings Model 

Perturbation theory breaks down in situation where we have a driving 
field of near resonance frequency. This is because the resonance causes large 
population transfers and the assumption made in perturbation theory that 
Cj (t) 1, no longer holds. Therefore we must take another approach. One 



such approach can be understood by noticing that for the case of near reso- 
nance, most of the population is transferred to the near resonant state, so that 
the other states may be neglected. The two most dominant states remain, and 
the resulting system of differential equations are much more tractable. 
The Jaynes-Cummings Hamiltonian is (see[l],[10]) 

+ aV) • (4.70) 

The interaction term is 

H int = h\ (a + a + a+a" 1 ) (4.71) 
and it induces the transitions, 

\e)\n) +*\g)\n + l) (4.72) 

or 

\e)\n-l) ^ \g)\n) . (4.73) 

The product states | e ) | n ) , \ g ) | n + 1 ) ,etc. , span what is known the bare 
basis, they are bare states of the Jaynes-Cummings model. For any given n, 
the dynamics of the system is confined to a 2D space of product states spanned 
by 

{|e)|n-l),|<7)|n-l),|e)|n),|<7)|n)}. (4.74) 



We can therefore define general product states for any n 



It follows that 



\ipin) = \e)\n) 



\^2n) = \g) \n + l) 



(4>ln\ \lp2n) = 



(4.75) 
(4.76) 

(4.77) 



The matrix representation of H 



Hij = (ip in \H\tp 



(4.78) 



In this basis becomes 



nu + \ fojjQ HXy/n + 1 
hXy/n + 1 (n + 1) uj — \hjjQ 



For any given n, the energy eigenvalues are 



E±(n)= [ n + - ) hto ± nn n (A) 



(4.79) 



(4.80) 



where A = (oj — oj), is the detuning parameter of the atomic transition fre- 
quency and the monochromatic field and 



n n (A) = [A 2 + 4A 2 (n+l)j 



1/2 



(4.81) 



is the damped Rabi oscillation frequency, which in the case of resonance, i.e, 
A = , becomes 

Q n (0) = 2Av / ^TT (4.82) 

The set of energy eigenstates form what are known as the dressed states, and 
these are a linear combination of the bare states which are 



\n, +) = cos ($ n /2) |V>m> + sin ($ n /2) |^ 2n ) 



(4.83) 



\n, -) = - sin ($„/2) \i/; ln ) + cos ($ n /2) |^ 2n ) 



(4.84) 



with $ n given by 



$ n = tan" 



! {2\y/n + T 



tan 



i (0) 



(4.85) 



Moreover, 



sin ($„/2) = 



cos ($„/2) 



v/2 



v/2 



(A) - A' 

a, (A) 

(A)+A 

Q n (A) 



1/2 



1/2 



(4.86) 



(4.87) 



The dressed states \n, ±) comprise the Jaynes-Cummings doublet. The Kl n (A) 
term splits the energies of the bare states] ipi n ) , \1p2n) , an effect known as the 
dynamic stark shift. In the case of exact resonance A = 0, the bare states 
become degenerate, but the dynamic stark shift splitting of the dressed states 
endures. In the exact resonance limit the dressed states can be represented in 



the basis of bare states as 



\n,+) = -j=(\e)\n) + \g)\n + l)) (4.88) 

\n,-) = ^={-\e)\n) + \g)\n + l)) (4.89) 
V2 

To obtain the dynamics in the dressed state basis, let us consider the case of 
an atom field system, for which the field is prepared in some superposition of 
initial states as 

foMO))=J]Cn|ra> (4.90) 

n 

and for which an atom, initially in the state |e) gets injected into the field. 
Therefore the initial state of the atom field system is 

1^/(0)) = hM0)> | e) (4.91) 

= Y J C nW)\e) = Y J °^n) ■ (4.92) 

n n 

Now, from (IV.2) and (IV.3), it follows that 

|Vm> = cos ($ n /2) \n, +) - sin ($ n /2) \n, -) (4.93) 
therefore the initial state of the atom-field system is 



IVV (0)} = °n ( cos ( $ -/ 2 ) K +) - sin ($„/2) \n, -)] (4.94) 

n 



One nice feature of the dressed states, is that they are stationary states of the 
Hamiltonian, and a consequence, the time evolution of H is 



= °n [cos ($ n /2) \n, +) e - lE+{n)t/h - sin ($ n /2) \n, -) e - iE ^ n)t/h ] (4.96) 



4.4 Jaynes-Cummings with large de-tuning 

In the forgoing we have been working with an "on resonance" approxima- 
tion, in which the detuning parameter vanishes i.e., when A = 0. A more 
interesting case is the one in which detuning exists to the extent that direct 
atomic transitions do not occur, but where dispersive interactions between a 
single and a cavity field do occur . 

The effective Hamiltonian in the case of large detuning is given by 




(4.95) 



n 



H e f f = kx [er+er_ + aW 3 ] 



(4.97) 



Where, 



X = A 2 /A 



(4.98) 



The transition operators are the projections 



°+ = |e> (g\ 



<J- = 



(4.99) 



Note that 

(7 +(7 _ = | e ) (e| (4.100) 
is the emission projector, and the inversion operator a 3 is 

a 3 = |e) (e| - |<7> (4.101) 

The transition and inversion operators obey the Pauli algebra 

[<7+,<7_] =<7 3 (4.102) 

[a 3 ,a ± ] = 2a ± (4.103) 

Suppose that the initial state of the atom-field system has the configuration 
of an atom in the ground state and the field in a number state, i.e. 

|^(0)> = b>|n> (4.104) 

The time evolved state becomes 

1^ (*)) = e -^//*M ^ (o)) = jxm |^ _ (4.105) 

Similarly, for the initial conditions 

|^(0)) = |e) \n) (4.106) 



We get 

1^ (f)) = e -^fft/h ^ ^ = e i X { n +i)t ^ ^ _ (4.107) 

and nothing happens except the production of unmeasurable phase factors. 

However, if the initial state is a coherent state of the field, that is in the 
case where 

(0)} = \g) \a) Coherent initial state (4.108) 

We obtain 

|^ {t)) = e-^ff^ \if> (0)) = \g) \ae ixt ) (4.109) 
Similarly, for the initial state 

|^(0)) = |e) |a) (4.110) 

we have 

\if, (t)) = e -*^//*/» |^ (o)) = e ~ lxt \e) \ae~ ixt ) (4.111) 

For either case of the initial coherent field state, the coherent state amplitude 
gets rotated in phase space by the angle 6 = xt- The direction of rotation 
depends on which initial states the atom is in. 

Let us now consider the case of an atom in an initial superposition of 
ground and excited states, which for simplicity, we assume takes the form 
of a balanced state with the form: 

l^atom) = 4= (\g) + e<* |e» phase (4.112) 



For an initial state 

|V>(0)) = |VWn) \<X) 



(4.113) 



We obtain 



llM*)> 



e -iH eff t/h ^ ^ 



1 



{\g)\ae lxt ) + e- t{xt -^\e) 



ixt )) (4.114) 



V2 



ae 



which is a much more interesting state, since now the atom and field are en- 
tangled. 

Taking xt = tt/2 , for which e %xt = i , and e~ %xt = —% , we have the entangled 
state 



which can be understood by analogy to the Schrodinger's cat paradox. With 
this analogy, our atomic states correspond to the radioactive atom in the para- 
dox, and the two phase-separated coherent field states play the role of Schrodinger's 
cat. Moreover the above entangled state corresponds the entangled state 

I "^atom-cat) = ~~j= (|atom not decayed) |cat alive) + |atom decayed) |cat dead)) 



Coherent states differing in phase by ir are maximally distinguishable, and 
there is effectively no overlap between the states, that is for |a| sufficiently 
large. Very large values of |a| are macroscopically distinguishable, while mod- 
erate values of | a |, mesoscopically distinguishable. 




(4.115) 



y/2 



(4.116) 



5 Application of CQED to Quantum Information 
Processing 



5.1 The Fabry-Perot Cavity 



We start this section with a brief overview the Fabry-Perot cavity. One of 
the most essential components of a Fabry-Perot cavity is a partially silvered 
mirror, which partially reflects and transmits incidents light E a and E b , which 
has the effect of producing output fields E a > and E v , which are related by the 
unitary transformation: 



E a > 




E^ 





y/T^R -y/R 





E a 




E b 



(5.1) 



where R is the reflectivity of the mirror. 

A Fabry-Perot (FP) cavity is made from two plane parallel mirrors of reac- 
tivities R\ and R 2 , incident upon which is light form outside the cavity E int . 
Inside the cavity light bounces back and forth between the two mirrors ac- 
quiring a phase shift e 1 ^ on each trip. The internal cavity field is 



5> 



V 1 — RE in 
1 + eVy/RJ^ 



(5.2) 



One of the most important things about the Fabry Perot cavity for purposes 
of CQED is the power in the internal cavity field mode as a function of of the 



power and frequency of the input field, 



P 



E 



2 



1 -Ri 



cav 



cav 



(5.3) 



in 



in 



1 + e^y/RJ^ 



Frequency selectivity arises because of constructive and destructive inter- 
ference between the cavity mode and the reflected light front. Another indis- 
pensable feature is that on resonance, the cavity field achieves a maximum 
which is approximately (1 — R)' 1 times the incident field . 

5.2 Quantum Computation 

Quantum information can be encoded with single photons in the dual rail 
representation 



Arbitrary unitary transformations can be applied to such quantum informa- 
tion using phase shifters, beam splitters, and nonlinear optical Kerr media 



The single photon representation of a qubit is attractive because it repre- 
sents the information saturation limit of the electromagnetic field and single 
photons, by today's standards can be generated relatively easily and more- 
over and most importantly, arbitrary qubit operations become possible, 
in general in the dual-rail representation. The difficult part in this ap- 
proach is making the photon - photon scattering amplitudes large enough for 
entanglement to occur. In an optical cavity this is commonly implemented with 



co|01) + ci|10) 



(5.4) 



[1]. 



optical nonlinear Kerr media. However, in reality even the best non-linear 
Kerr media are weak, and are unable to provide a cross phase modulation of 
180° between single photon states. It is estimated , that even in the best cases, 
approximately 50 photons would have to be absorbed for each 180° cross phase 
modulated photon. 

Despite its drawbacks, the optical quantum computer does provide us with 
some insight into the architecture and design of a quantum computer. As- 
suming we had sufficiently good components available, we could construct an 
optical quantum computer, which will be almost entirely comprised of opti- 
cal interferometers. The quantum information is encoded in both the photon 
number states and the photon phase. Interferometers perform the function 
of switching between the two representations. Stability however, becomes a 
major issue, and if a massive representation of a qubit is chosen, then sta- 
ble interferometers would be a challenge to construct because of the relatively 
short scale of the de Broglie wavelengths of the qubits. 
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